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Abstract. We prove higher-order and a Gevrey class (spatial analytic) reg- 
ularity of solutions to the Euler-Voigt inviscid a-regularization of the three- 
dimensional Euler equations of ideal incompressible fluids. Moreover, we es- 
tablish the convergence of strong solutions of the Euler-Voigt model to the 
corresponding solution of the three-dimensional Euler equations for inviscid 
flow on the interval of existence of the latter. Furthermore, we derive a cri- 
terion for finite-time blow-up of the Euler equations based on this inviscid 
rcgularization. The coupling of a magnetic field to the Euler-Voigt model is 
introduced to form an inviscid regularization of the inviscid irresistive magneto- 
hydrodynamic (MHD) system. Global regularity of the regularized MHD sys- 
tem is also established. 



This work is dedicated to Professor Peter Kloeden 
on the occasion of his 60th birthday. 

1. Introduction 

One of the outstanding problems in mathematical physics is to find an accurate, 
practical description of turbulent flows. This problem is not only out of reach for 
current mathematical tools, but direct numerical simulation of detailed turbulent 
flows has proven to be computationally prohibitive. This is due to the current 
inability to resolve the wide range of the underlying spatial and temporal scales, 
even by using the most powerful state-of-the-art computers. For the case of an 
incompressible, viscous, homogeneous fluid in a domain fi C I 3 , the governing 
equations are widely thought to be given by the Navier-Stokes equations 

(1.1a) d t v + (v ■ V)v + Vp = vAv + f infix [0,T), 

(1.1b) V-u = infix [0,T), 

(1.1c) w(0) = v m in fi, 

(l.ld) periodic, or v = on dfl. 

where v(x,t) = (vi(x,t),V2(x,t),V3(x,t)) denotes the velocity field of the fluid at 
the point x = (x\, X2, X3) G fi and the time t € [0, T), T > 0; p(x,t) denotes 
the pressure; f(x,t) = (fi(x, t), f2[x, t), fa{x, t)) is the body (external) forcing on 
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the fluid and v > is the kinematic viscosity. The Dirichlct boundary conditions 
v\fm = correspond to the physically relevant no-slip case, while the equations un- 
der periodic boundary conditions are mathematically simpler, and preserve many 
(though not all) of the structures present in the Dirichlet case. In the inviscid case, 
which is the focus of this paper, the equations are known as the Euler equations. 
These equations are identical to (1.1), except that v = and the boundary condi- 
tions (1-ld) either remain periodic, or are replaced by the Neumann no penetration 
boundary conditions v ■ n = on dfl, where n is the outward point normal vector 
O. Despite much work on the three-dimensional Euler equations over the last two 
and a half centuries, many basic questions remain unanswered. For recent surveys 
of the known results about the Euler equations, see, e.g., [6,16,42]. 

As a result of the unresolved difficulties mentioned above, researchers have fo- 
cused on obtaining reliable, computable models for the large-scale behavior, by 
considering, for example, spatial or temporal averages of the Navier-Stokes equa- 
tions. However, because of the nonlinear nature of (1.1), it is not possible to obtain 
an exact closed analytical system for the large-scale motion. Thus, one must re- 
sort to modeling the large-scale motion of the dynamics of (1.1) for turbulent flows 
without having to compute their full-scale dynamics. In many systems, this is often 
achieved by filtering, or taking a localized spatial averaging, of the equations for a 
quantity tp to get equations for some filtered quantity Tp. However, if we apply this 
method to the non-linear equations (1.1) (say, under periodic boundary conditions), 
we arrive at 

(1.2a) d t v + V • (vWv) + Vp = v/\v + J, 

(1.2b) V-U=0, 

where we have used the fact that (v ■ V)w = V • (v ® v), due to the incomprcssibility. 
Thus, since v <8> v ^ v®v for any useful definition of the filtering v, (1.2) is a system 
of equations in both v and v (or, equivalcntly, for the mean v and the fluctuation 
v — v), rather than just v alone. The fact that we now have more unknowns than 
equations is known as the closure problem of turbulence. By adding and subtracting 
V- (v<& v) in (1.2a), we realize that we must deal with the Reynolds stress tensor, 

TZ(v, v) := v ® v — v (g> v. 

Much effort by researchers modeling turbulence has gone into finding useful ap- 
proximations to TZ(v, v) in terms of v alone. 

The approximation lZ(v, v) « v ®v — v ®v was introduced and studied by 
Bardina in [3]. Later, in [36], a simpler approximation was considered, namely 
TZ(v, v) k v (giv — v ®v, where in particular the Helmholtz filtering operator Tp := 
(I — a 2 /\)~ 1 (p was used, where a > is a length scale that represents the width 
of the spatial filter. (Note that here, the inverse is taken in the context of mean- 
zero functions with periodic boundary conditions). This choice of filtering has 
proven to be very effective in turbulence modeling (see, e.g., [11-15,24,27]), in 
particular, it is a key ingredient to successes of the a-models of turbulence, which 
were introduced as a contribution to the efforts of modeling large-scale motion. 
These analytical models have shown much promise in recent years, notably in the 
benchmark cases of turbulent flows in pipes and channels [11-15,24,27]. Namely, 
analytical solutions to the filtered steady-state equations of the a-models were found 
to successfully match time-averaged experimental data of turbulent flows for a wide 
range of large Reynolds numbers (see, e.g., [11-15]). Using the Helmholtz filtering 
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and the approximation lZ(v,v) « v(£iv — v ®v in (1.2) results in the following 
model, known as the simplified Bardina model: 

(1.3a) (/ - a 2 A)d t u + (u ■ V)w + Vp = vA{I - a 2 A)u + /, 

(1.3b) V-u = 0. 

where u := v = (I — a 2 A)~ 1 v. Note that in the case where a = 0, (1.3), as a 
PDE, coincides formally with equations (1.1) under periodic boundary conditions, 
a fact which is true for all other a-models (see, e.g., [12-15,24,27]). Model (1.3) 
was investigated in [11,36]. In [11], the authors proved the global existence and 
uniqueness of solutions in both the viscous (y > 0) and inviscid (y = 0) cases. In 
particular, the authors prove in [11] that for initial data u m £ H 1 ^), system (1.4), 
below, has a unique solution u G C 1 ((— oo, oo), i/ 1 (0)). This result is of particular 
interest in the inviscid case since, as of yet, (1.4) is the only a-modcl for which 
global regularity in the inviscid case has been proven. Furthermore, it was noted 
in [11] that formally setting v = in (1.3) amounts to simply adding the term 
—a 2 dtAu to the left-hand side of (1.1a) (with v = 0), yielding 

(1.4a) -a 2 d t Au + d t u + (u ■ V)it + Vp = /, 

(1.4b) V • u = 0, 

(1.4c) u(0) = u m , 

which we call the Eulcr-Voigt equations. The boundary conditions are taken to be 
periodic, and the domain is the periodic unit torus, f2 = T 3 := [0, l] 3 . Furthermore, 
we impose the mean-zero condition, 

(1.5) / u m dx= [ fdx = 0, 

which implies that J T3 u dx = 0. 

Remarkably, if one reintroduces the viscous term vAu to the right-hand side of 
(1.4a), the resulting equations happen to coincide with equations governing certain 
visco-elastic fluids known as Kelvin- Voigt fluids, which were first introduced and 
studied by A. P. Oskolkov [46,47]. These equations are known as the Navier-Stokes- 
Voigt equations, which is our reason for calling (1.4) the Euler- Voigt equations. 
They were proposed in [11] as a regularization for cither the Navier-Stokes (for v > 
0) or Euler (for v = 0) equations, for small values of the regularization parameter 
a. 

In the presence of a physical boundary, and under the assumption of the no-slip 
boundary conditions u\dQ = 0, the Navier-Stokes- Voigt equations (i.e., system (1.4) 
with vAu added to the right-hand side of (1.4a)), as an a-model regularization for 
the three-dimensional Navier-Stokes equations, have a very attractive advantage 
over other a-models and subgrid-scale models in that one does not need to impose 
any additional artificial (i.e. non-physical) boundary conditions to prove the global 
existence and uniqueness of strong solutions, as it has been pointed out in [11]. 
It is worth mentioning that the long-term dynamics and estimates for the global 
attractor of the three-dimensional Navier-Stokes- Voigt model was studied in [30]. 
Moreover, it was shown recently in [50] that the statistical solutions (i.e., invari- 
ant probability measures) of the three-dimensional Navier-Stokes- Voigt equations 
converge, in a suitable sense, to a corresponding statistical solution (invariant prob- 
ability measure) of the three-dimensional Navier-Stokes equations. Furthermore, in 
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the context of numerical computations, the Navier-Stokes-Voigt system appears to 
have less stiffness than the Navier-Stokes system (see, e.g., [20,37]). 

In [37], the statistical properties of the Navier-Stokes-Voigt model were inves- 
tigated numerically in the context of the Sabra shell phcnomcnological model of 
turbulence and were compared with the corresponding Navier-Stokes shell model. 
In particular, it was observed that for values of the regularization parameter a 
smaller than the Kolmogorov dissipation length scale, the structure functions of 
both models obey the same power laws in the inertial range. For values of a larger 
than this scale, two distinct regions associated with the inertial range of the en- 
ergy spectrum for the Navier-Stokes-Voigt model were observed in [37]; namely, a 
region of low wave numbers obeying the Kolmogorov fc 2 / 3 power law, and a region 
of higher wave numbers, where energy condensates. 

Due to its simplicity, the Voigt a-rcgularization is also well-suited to being ap- 
plied to other hydrodynamic models, such as the two-dimensional surface quasi- 
geostrophic equations, demonstrated in [33], and the three-dimensional magncto- 
hydrodynamic (MHD) equations, which we demonstrate in this contribution. It is 
also worth mentioning that in the case of the inviscid Burgers equation, ut+uu x = 0, 
this type of regularization leads to —a 2 u xxt + u t -\-uu x = 0, which is the well-known 
Bcnjamin-Bona-Mahony equation of water waves [8]. 

In the present paper, we focus on the Euler- Voigt equations subject to periodic 
boundary conditions. Our results are organized as follows. In Section 2, well-known 
results are stated and standard notation is recalled. In Section 3, we extend the 
results of [11] to prove higher-order regularity for (1.4) under periodic boundary 
conditions. In Section 4, we prove that the solutions enjoy spatial analyticity for 
u m analytic. In Section 5, we prove that solutions to the Euler- Voigt equations 
converge, in some sense, to sufficiently regular solutions of the Euler equations, 
as a — > 0, on any closed interval of time where the corresponding solutions of 
the three-dimensional Euler equations exist. Moreover, we provide a criterion for 
blow-up of the Euler equations, based on the Voigt a-regularization. In Section 
6, we demonstrate the applicability of the Voigt a-regularization to other fluid 
models by considering a similar a-regularization of hydrodynamic models involving 
magnetism. In particular, we prove the global existence of solutions for the inviscid 
irresistive MHD- Voigt regularization. The theorems given in this paper also hold 
in K 3 (see, e.g., [45]), but for simplicity we work in T 3 . Our methods and results 
readily apply to the viscous (i.e., Navier-Stokes-Voigt) case as well. For simplicity, 
we set the forcing / equal to zero in all models considered, although our results 
hold with little additional effort, given a suitably smooth forcing term satisfying 
(1.5). 

2. Preliminaries 

In this section, we introduce some preliminary material and notations which are 
commonly used in the mathematical study of fluids, in particular in the study of 
the Navier-Stokes equations (NSE). For a more detailed discussion of these topics, 
we refer to [18,25,54,55]. 

Let T be the set of all vector-valued trigonometric polynomials with periodic 
domain T 3 := [0, l] 3 . We define the space of 'test' functions to be 

V := j</? S T : V • ip = and J ip(x) dx = j . 
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We denote by LP and H m the usual Lebesgue and Sobolev spaces over T 3 , and 
define H and V to be the closures of V in L 2 and H 1 respectively. We define the 
inner products on H and V respectively by 

(u, v) = S^ [ u t Vidx and ((u,v)) = / ^-^-dx, 
^■1 ■ ij^i^ 3 dx i dx 3 

and the associated norms \u\ = (u, it) 1 / 2 , ||u|| = {{u^u)) 1 / 2 . Note that ((-,-)) is a 
norm due to the Poincare inequality, (2.2), below. We denote by V' the dual space 
of V. The action of V on V is denoted by (•, •) = (•, -) v ,. Note that we have the 
continuous embeddings 

(2.1) V^H^V'. 

Moreover, by the Rellich-Kondrachov compactness theorem (see, e.g., [1,21]), these 
embeddings are compact for bounded domains. 

Let X be a Banach space. We denote by L p ((a, b),X) the space of Bochner 
measurable functions t i— > w(t), where w(t) £ X for a.e. t £ (a, b), such that 
the integral ||w(*)llx dt is finite (see, e.g., [1]). A similar convention is used for 
C fc ((a, b),X). Abusing notation slightly, we write w(-) for the map t n> w(t). In 
the same vein, we often write the vector-valued function w(-,t) as w(t) when w is 
a function of x and t. 

We denote by P a \ L 2 — y H the Leray-Hclmholtz projection operator (i.e., the 
orthogonal projection onto solenoidal, i.e., divergence- free, vector spaces), and de- 
fine the Stokes operator A := —P a A with domain V (A) := H 2 fl V. In our case 
of periodic boundary conditions, it is known that A = —A (see, e.g., [18,54]). 
A^ 1 : H — >• H is a positive-definite, self-adjoint, compact operator, and therefore 
has an orthonormal basis of cigenfunctions ip^ corresponding to a non-increasing 
sequence of eigenvalues (see, e.g., [18,54]). We observe that (I + a 2 A) -1 is a well- 
defined bounded operator and that (I + a 2 A) and (/ + a 2 A)" 1 arc self-adjoint. 
Furthermore, due to the periodic boundary conditions, partial derivatives of any 
order commute with (/ + a 2 A) and (/ + a 2 A) -1 . We label the eigenvalues of 
A so that < Ai < A2 < A3 < • ■ • . Notice that in the case of periodic boundary 
conditions in the torus T 3 we have Ai = (27r)~ 2 . Furthermore, for all w 6 V, we 
have the Poincare inequality 

(2.2) HHU=(T3) < A i ; 1/2 ||Vw|| i 2 (T 3 ) . 
Due to (2.2), for w £ T>(A), we have the norm equivalences 

(2.3) \Aw\ = |H| ff 2 and || Vu;|| L 2 (T 3 ) ^ \\w\\. 
It will be convenient to use the standard notation 

(2.4) B(wi,w 2 ) :=P a ((wi ■ V)w 2 ) 

for w\,W2 £ V. We list several important properties of B which can be found for 
example in [18,25,54,55]. 

Lemma 2.1. The operator B defined in (2.4) is a bilinear form which can be 
extended as a continuous map B : V x V — > V . Furthermore, for w\, W2, W3 £ V, 



(2.5) 



(B(w 1 ,W2),w 3 ) v , = - (B(wi,W3),u} 2 ) v , , 
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and therefore 

(2.6) (B(w 1 ,w 2 ),w 2 ) v , = 0. 

Here and below, C, c,-, etc. denote generic constants which may change from line 
to line. C a ,C' a7 C(- ■ ■ ), etc. denote generic constants which depend only upon the 
indicated parameters. 

Next, we recall Agmon's inequality (see, e.g., [1,2,18]). For w £ D(A) we have 

(2.7) ||HU~(n) <tf|M| 1/2 PH Va . 

Finally, we note a result of deRham [55,56], which states that if g is a locally 
integrable function (or more generally, a distribution), we have 

(2.8) g = Vp for some distribution p iff (g, v) = Vv £ V. 

Using this theorem, it can be shown (see, e.g., [18,55]) that (1.1) with v = 0, / = 0, 
and periodic boundary conditions, is equivalent to the functional equation 

du 

(2.9) — + B(u,u) = 0, 
and that (1.4) is equivalent to 

(2.10) (I + a 2 A)^ + B(u,u) = 0. 

at 

where the last equality is understood to hold in the sense of V . Finally, we define 
the notion of a solution to (2.10) which was given in [11]. 

Definition 2.2. Let u in £ V and consider a time interval (T 1} T 2 ) with Ti < < T 2 . 
A function u 6 C 1 ((Ti, T 2 ), V) is said to be a solution to (2.10) if it satisfies (2.10) 
in the sense of V and furthermore u(0) = u m . 



3. H m Regularity 

In this section we show that solutions to (1.4), or cquivalently (2.10), are globally 
well-posed and have as much smoothness as the initial data. However, we do 
not expect the solutions to instantaneously gain additional smoothness, even in 
the viscous (Navier-Stokes-Voigt) case. The reason for this is that adding the 
term —a 2 d t Au to the Navier-Stokes equations destroys their parabolic structure 
(indeed, the Navier-Stokes-Voigt equations are well-posed backwards in time [11], a 
feature which is not present in any parabolic equation). Instead, the Navier-Stokes- 
Voigt system behaves like a damped hyperbolic (pseudo-parabolic) system [30]. 
Despite this, the Navier-Stokes-Voigt equations (with analytic forcing) possess a 
finite dimensional global attractor comprised of analytic functions. For details on 
these matters, we refer to [28-30]. With these considerations in mind, we now state 
the following theorem. 

Theorem 3.1. Let u™ £ H m (T 3 ) n V, for m > 1. Then there exists a unique 
solution u of (2.10) (with f = 0, for simplicity) with u £ C 1 ((— oo, oc), iJ m (T 3 ) n 
V). Moreover, 

IK*)||h-(t3) < C(a, \\u in \\ Hm )(l + \t\Y* m \ 
for all t £ (— oo, oo), where p(l) = 1, p(2) = 2 and p(m) = 5 (§)"' 3 — 1 for m > 3. 
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We give two proofs of Theorem (3.1). Proof 1 is based on the contraction method, 
which demonstrates the ODE nature of (1.4). Proof 2 uses the Galcrkin approxima- 
tion procedure, which is essential for proving Theorem 4.2 below, and furthermore 
is of interest from the point of view of numerical analysis (see, e.g., [18,25,43,54,55] 
for detailed discussions of this method). 

Proof 1. The existence and uniqueness of solutions in the sense of Definition 2.2 
have been established in [11]. We prove only the higher-order regularity. Applying 
(7 + aA)^ 1 to (2.10) (the inverse taken with respect to the periodic boundary 
conditions), we obtain 

(3.1) ^ = {I + a 2 A)~ 1 B{u 1 u) :=N(u). 
dt 

We claim that N : H m n V — > H m n V is locally Lipschitz continuous for each 
m > 1. The case m = 1 was proven in [11]. We prove the remaining cases to > 2. 
Let ui,U2 € H m n V be arbitrary and write 5u = u 2 — u\. We have 

\\N{ Ul ) - N(u 2 )\\ H «, 
= SU P (9 (I + a 2 A)- 1 (B(u 1 ,u 1 )-B{u2,u 2 )) > d ( p) 

0<\/3\<m IMI» m=1 

= SU P (^(S(fe, Ml ) + B(i t2 , ( 5M)),a' 3 (/ + a 2 A)-V) 

0<|/3|<m \\V\\ H ™=1 

= sup V r) (B{d' i 5u,d p -^u 1 ) + B{d'-'u 2 ,d f} -' i 5u),d l3 {I + a 2 A)- 1 ip) . 
IMI*»= 1 o<feA 7 >' 

0<7</3 

In the case 7 = 0, we have for |/3| < m, 

{B{5u,d u 1 ),d p {I + a 2 A)- 1 ip) = - (B(Su, d@(I + a 2 A)^ 1 Lp) , d 13 u\) 

< ||^|| LO c (T 3 ) ||V^(/ + « 2 yl)-V!l^(T3)||^ U i|| i2(T 3 ) 

< Ca||<MlH 2 (T3)IMIff"— i(T3)||ui||i?"»(T 3 )- 

Similarly, 

(B(u a ,d fi Su),8 p (I + o/'A)- 1 (p) ^Calluallfla^H^llfl^-i^llH^CT"). 
In the case \j\ > 0, we have for |/3| < to, 7 < /3, 

(5(<9 7 <5u, 0"(J + a 2 A)-V) 

< ||^fe|| i2(T 3 ) ||V^-T U i|| i2(T 3 ) ||^(7 + a 2 A)-V|U«>(T3) 

< C a ||5u||ijm( T 3)||ui ||ifm( T 3) ||^||ffm( T 3). 

Similarly, 

(Bfff'uj.^ttj/g + ^-V) <C a ||u a || fl3(X 3 ) ||HlH»ci a )llvllfl«(i»)- 
Combining the above estimates, we have for m > 2, 

(3.2) ||J\T(ui) - JVMHirm < C a , TO (||u 1 ||h-»cp«) + ll«a||n*(T»)) HHIj?»CF«)- 

Choose i? such that ||wi||_fr">(T 3 ), ||w2||_f/™(T 3 ) < 7?. Then (3.2) shows that N is 
locally Lipschitz in H m (T 3 )DV , with Lipscitz constant 2RC a , m in the ball of radius 
R centered at the origin. Let u m € H m (T 3 ). By the Picard-Lindelof Theorem, there 
exists a time T > such that (3.1) has a solution in C 1 ([-T, T], F m (T 3 ) n V). Let 
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[0,T^) be the maximal positive interval of existence. (For the sake of clarity, 
we work only in the positive time case, t > 0; however, the same proof holds 
for t < 0.) All of the work below takes place on the interval [0,T^). If T* n < 

00, then limsup t ^ r « ||it||jy m = oo, but as we will show below, u is bounded in 
C([0, T,*), H m (T 3 ) n V), and thus we will have T,* n = oo. The remainder of the 
proof is divided into four steps. The cases m — 1, 2, 3 are treated sequentially. The 
general case is then shown inductively. 

1. For the case m = 1, we take the inner product of (1.4a) with it, integrate by 
parts, and employ (2.8) to arrive at 

f« 2 ||Vw|||2( T 31 + ]]u|||2 fT 3^ =0, 



~2dt V 11 llL2 ( T3 ) ll""£ 2 (T 3 ) y 
Integrating this inequality gives 

(3.3) a 2 \\Vu(t)\\ 2 L2{T3) + H*)|li3(T3) = a 2 ||Vuo||i2 (T3) + IKHl^t 3 ), 

which implies || 'u(i) || i^ 1 (T 3 ) < Cq = C(a, ||u m ||), and thus Tj* = oo. Note that 

(3.3) was essentially obtained in [11] to show global existence of solutions to (1.4) 
(equivalently (2.10)). 

2. For the case m = 2, we take the inner product with -Ah and integrate by parts 
to obtain 

(3.4) ~ (a 2 || Au\\ 2 L2{T3) + || V U || 2 2(T3) ) - ((u ■ V)u, Au) , 

where we have used (2.8). We estimate the right-hand side using the Holder's 
inequality and Agmon's inequality. 



| ((u • V)u, Am) I < / M|Vu||Au| dx 

< lk||L~(T 3 )l|Vu|| L 2 (T3) ||Au|| i 2 (T3) 

< C||W||^1 2 (T3) ||W||H2 (T 3) II Vu||i2 (T 3) || Au|| ia(T 3) 

— C|MliJi(T 3 ) ll U ll_f/ 2 (T 3 ) 

< C Q ,||u|| i3 - 2 ('J3)) 

thanks to (3.3). Using this estimate in (3.4) and employing (2.3) gives 

d ( \ / \ 3 / 4 

(3.5) - [a 2 \\Au\\l 2(T3) + HVtilliafpjJ < C a (a 2 ||A w || 2 2(T3) + ||V W || 2 2(T3) J . 

Thus, using Gronwall's inequality and the norm equivalence (2.3), we have the 
algebraic growth rate in time 

(3.6) IKi)||ff 2(T 3) < C a (||tf"||^ 2 (T 3) +t)\ for all t e [0,T\. 
Thus, T 2 * = 00 • 

3. For the case m = 3, we take the inner product with A 2 u, integrate by parts and 
use (2.8) to obtain 

(3-7) ™ (a 2 ||VAu|| 2 2(T3) + ||Au||i a(x8) ) =-((«• V)u, A 2 u) . 
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For the right-hand side, we integrate by parts twice to obtain 



/ [(u • V)u] • A 2 udx 

JT 3 

A((u ■ V)u)] ■ Audx 
(Au ■ V)u + (u • V) Aw + 2 diii ■ V<% 

[ [{Au ■ V)u] • Ait dx + 2 V / (<%u • V<9,u) • Au 

JT 3 ~{ JT 3 



Au dx 



WAi 



Il 2 (t 3 ) 



|Vu| 



L oo( T 3 



2||Vu|| L oo (T 3 ) ||VVw|| i 2 (T 3)||Au|| i 2 (T 3 ) 



<C||n||l, 2(T 3 ) ||Vu|Uo 0(T 3 ) 



2 (T 

<C|| U || 2 ff 2 (T 3jV U 



,1/2 

Lf^Cr 3 ) 



IVii 



1/2 

Lf/ 2 (T 3 ) 



— C|l U ll_f/2( T 3) ll U ll#3( T 3) 



on [0, T]. For the third equality above, we have employed the fact that 
J T3 ((u ■ V)Ait) • Audx = 0, which follows from (2.6) since V • u = 0. The first 
inequality is due to Holder's inequality, and the third is due to Agmon's inequality 
(2.7). Using the above estimate and (3.6) in (3.7) gives 



(3.8) j t (a 2 ||VAw||2 2(T3) + ||A 



U\\~ L 2(J3) I < 



C' a (c a t+\\u m \\\ 



I J/ 3 (T 3 ) • 



For the last factor on the right-hand side of (3.8), one can use the norm equivalence 
(2.3) to show that 



(3.9) 



l u ll// 3 (T 3 ) — 



C a (a a ||VAu||i a(T , ) + ||Au||| ! i (T , ) ) 



Combining (3.8) and (3.9), we may apply Grdnwall's inequality and again the norm 
equivalence (2.3) to find the algebraic growth rate in time 



IM| ff 3 (T 3 ) < IKl^T 3 ) + C '^ Cat + H Um |l^(T 3 )) £ 



2/3 



(3.10) 



<C a 



\u m \\ H3iJ3) + (t + \\u m \\ H2iT3) y 



Thus T 3 * = oo. 



4. For the remaining cases m > 4, we work by induction on m. That is, we show 
that if the iJ m_1 norm is bounded, then so is the H m norm. We follow closely 
the techniques used in [43] (see also [42,53]). Let (3 be a multi-index such that 
|/3 1 < m — 1. We apply d& to both sides of (1.4a), take the inner product of (1.4a) 
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with d^u and use (2.8) to obtain 



I|( a 2 ||v ^ uj| 2 2(T3) + || ^ u|| 2 2(T3 

= - V)u),d l3 u) . 



(3.11) =- V (P) [ \{cTiu-V)d^u]-d udx 

since in the case 7 = we have J J3 [(u ■ V) d^u] ■ d^u dx = 0, again by the divergence 
free condition (2.6). We estimate the integrals in (3.11), analyzing the cases I7I = 1, 
I7I = 2 and I7I > 3 separately. 
For I7I = 1, 7 < /3, we have 

T 3 

< CIMIimi+slMl^i =C||u||h»||«||^,. 

For I7I — 2, 7 < f3, we have 

/ p 7 u- v)a /3 ^ 7 u] -d p udx < ||a 7 u|| L 6 (T 3 ) ||va /3 - 7 u|| L3(T 3)||^u|| i 2 (T 3 ) 

<C\\d^u\\ m{T s ) \\Wd^u\\ HHT3) \\d^u\\ L 2 {T3) 

< C||lt||_ fJ -3 (T 3 ) ||u|| ff |3|-l( T 3)||M|| // |jJ| (T 3 ) 



< 



T 3 



For I7I > 3, 7 < /?, we have 

[(d-y u ■ v) d^u] ■ d^udx < Wd-'uW^^wvd^uW^^Wd^uW^^ 

t 

< C|l«llHlTl(T 3 )l|V9 /3 " 7 u|| ff 2 (T 3 ) ||a' 3 u|| i 2 (T 3 ) 

< C||«llfflTl(T 3 )IMIffl/3|(T 3 )- 

These estimates imply that the right-hand side of (3.11) is bounded above by 
C(|M|ff 3 (T 3 ) + |M|ffifi(T 3 ))IMlHi3i(T 3 ) - C \\ u \\h™-i(t3)- Using this with (3.11) and 
summing over all j3 with (3 \ < m — 1 gives 

(3.12) j t (a 2 ||Vu||^ m _ 1(T3) + \\u\\ 2 Hm - 1[T3) ) < C\\uf Hm - 1(J3y 

Integrating this inequality and using the norm equivalence (2.3), we find the fol- 
lowing recursive relationship for m > 4: 

(3.13) \\u(t)\\ Hmm < C a (VlW 3 ) + / IK s )llH--UT 3 ) ds ) • 

Thus = co. Repeatedly iterating estimate (3.13) and eventually using (3.10), 
we see that ||u(i)||£f m ( T3 ) grows at most algebraically in t, that is 

(3.14) \\u(t)\\ Hm(j3) < C(a, |k n ||iP*(T3))(l+t) p(m) , for all t 6 [O.oo). 
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for some algebraic growth rate p(m). Using induction and (3.13), we have p(m) — 
(3p(m — 1) + 1) /2 for m > 4. From (3.10), we havep(3) = 4. Solving this difference 
equation, we find p{m) = 5 (3/2)'™ — 1, as claimed. □ 

Proof 2. For each N G N, let Pn denote the orthogonal projection in H given by 
Pn ■ H — > span {<Pk}\ k \ <N ■= Hn (fk are defined in Section 2). The Galerkin 
approximation to (2.10) at level N is given by the following system: 

(3.15a) (/ + a 2 A)^-u N = -P N B(u N , u N ), 

at 

(3.15b) u N (0) =P N u in . 

Applying the operator (/ + a 2 ^4) _1 to (3.15a), we see that this is an ODE in 
the finite-dimensional space Hn with quadratic non-linearity, and therefore has 
a solution un G C 1 ((— 7jv, Tn), Hn) for some Tn > 0. Thus, taking the inner 
product of (3.15a) with un and integrating by parts is justified, and we arrive at 

(3.16) a 2 ||Vuiv(t)||l2 (T3) + ||uAr(*)||i2 (T 3 ) < af 2 ||Vuo|| 2 2 (X 3) + lluolli^xa)- 

This implies that Tn = oo for all N > 1. Let T > be fixed but arbitrary. For 
simplicity, we work on the interval [0,T], but the same proof holds on [— T, 0]. 
Following ideas similar to the case of the 3D Eulcr equations (see, e.g., [43]), we 
show that {un} N£N is a Cauchy sequence in C([0,T], V). Let u m G V, and let un 
and um be solutions of (3.15) for N, M G N, N < M respectively, and consider 
uu M := un — um G Hm- We will need the following inequalities, which follow easily 
from Parseval's identity. 

||ujv||#2( T 3) < CN\\u N \\m(T3) and - Pn)um\\l^) < CN~ X \\u M \\m(j^)- 

Subtracting the equations (3.15a) for un and Um, applying / + a 2 A to both sides, 
taking the inner product with w M and using (2.5) and (2.6), we find 

\db (ll u; ^l^ 2 ( T3 ) + a2 H Vw Mlli 2 (T3)) 
= (B(un, un), (I - Pn)um) - (B(w M , un),w m ) 

< C|kjv||ffi 2 (T3)ll"Jv||ff2 2 (T 3)||(-f - Pn)um\\l*{t*) + C\\wm\\w(tz)\\ u n\\h1(t*) 

< C||ltj V ||^ / i 2 (T 3 ) A rl/2 ||MAr||^ / i 2 (T 3 ) A r_1 ||wAf||i/i(T3) + C\\ W M \ \ jji (T 3) || U N \\ H i (T3) 

< C a N-^ 2 + C a (\\w M \\h {r3) +« 2 ||V<||i 2(T3) ) , 

since ||itjv ||^f 1 (T 3 ) < Co uniformly in N. Grdnwall's inequality yields 

\\w M (t)\\h { ^)+^\\^M(t)\\h(^) 

< (lK(0)lli 2(T3) + a 2 \\V(w M (0))\\h m ) e C °' + C^-^e -* 1). 

Since u m G V, it follows that \\w M (0)\\ 2 L2{T3) +a 2 \\V{w M {0))\\ 2 L2{T3] as M, N 
oo. Therefore, the above inequality shows that {un}^^ is a Cauchy sequence and 
hence converges to an element u G C([0,T],V). Choose v £ V arbitrarily. Since 
u« G C([0, T],V), taking the inner product of (3.15a) with v and integrating in 
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time is justified, so we have 

(u N (t), v) + a 2 {{u N (t), »)) - (u N (0), v) - a 2 ((u N (0), v)) 

(3.17) =-/ (B(u N (s),u N {s)),P N v)ds. 

Jo 

Since un — > u strongly in C([0,T] f V), we have (uN(t'),v) — > (u(t'),v), and 
((u N (t'),v)) -> ((u(t'),v)) for all t' G [0,T]. Furthermore, we note that 

(B(un(s),un(s)), Pnv) ds — / (B(u(s),u(s)),v) ds 

Jo 

t pt 
(B(u N (s),u N (s)) 1 P N v - v) ds + / (B(un(s) — u(s),u N (s)),v) ds 

Jo 

t 

(B(u(s), un{s) — u(s)), v) ds. 

Using this and the facts that H^jvll jyi(T 3 ) is bounded independantly of N, that 
un — > u strongly in C([0,T], V), and that \\PjfV — v Wh 1 ^ 3 ) —> as TV — > oo, we find 
that Jq(B(un(s),un(s)), Pnv) ds — » f Q (B(u(s),u(s),v) ds as A~ — > oo. Passing to 
the limit in (3.17), we have 

(ii(t),w)+a 2 ((u(t),w))-( M m ,w)-a 2 ((u m ,w)) = - f (B(u(s),u(s)),v) ds 

Jo 

for all v G V. A simple argument using the density of V in V shows that the equality 
holds for all v G V. Since u G C([0,T],V) we have Au, B(u,u) G C([0, T], V). 
Therefore the above equality shows that (2.10) holds in the sense of V and that 
u(0) = u ln . Thus, we have proven the existence of a solution to (2.10). Thanks 
to (2.10) and the fact that {I + a 2 A)- 1 B{u,u) G C([0,T],V), we have f t u G 
C([0,T],V) so that u G C 1 ([0,T],V). To show the uniqueness and continuous 
dependence on initial data, let Ui,u 2 G C([0,T],V) are two solutions of (2.10) 
initial data u™,it 2 n G ^ respectively, and we write Su := u 2 — u\. Subtracting, we 
have 

(3.18) ^-8u = (I + a 2 A)- 1 (-B(u 2 ,Su) - B(5u, Ul )). 
dt 

Since jj^Su G C([0, T], V), we may justifiably take the inner product of (3.18) with 
5u and move the time derivative outside the inner product. Using (2.5) and the 
fact that {I + a 2 A) -1 is self-adjoint, we have 

^4ll*«lli»(T») = {B(u 2 ,(I + a 2 A)- 1 5u),6u) - (B(Su,ux),(I + a 2 A)- x 5u) 

< ||«2||l«(T3)||V(J+ Q! 2 A) _1 (5u|| L 3( T 3 ) ||(5M|| i 2( T 3 ) 

+ 11^1^2(^)11^7^11^2(^)11(7 + a 2 AY 1 5u\\ L ^ [T 3 ) 

< C a (\\u 2 \\hi(t 3 ) + ||wi||jji( T 3)) ||<Hli*(T3) ^ C a\\Su\\l 2( j 3 y 
since |ui Hi^i (t 3 ) 7 ||w2||ffi(T 3 ) < C a . By Gronwall's inequality, we see that 

\\u 2 {i) - Mi(t)||i2( T 3) < \\u 2 n - ul"|| 2 j2 ( T 3)e c '" t . In particular, if u" % = u\ n , then 
u 2 (t) = ui(t) for all t G [0,T]. 

Finally, we show the higher-order regularity. Suppose that u m G H m (T 3 ) n V. 
The process leading to (3.6), (3.10), and (3.14) is justified for un, since un G 
C([0, T], V). Thus, u N is bounded uniformly with respect to N in C([0, T], iJ m (T 3 )n 
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V). By the Banach-Alaoglu theorem, we see that un has a subsequence con- 
verging in the weak-* topology to an element v G L°°([0, T], H m (T 3 ) n V), and 
furthermore that the estimates (3.6), (3.10), and (3.14) hold for a.e. t G [0,T]. 
On the other hand, we already know that un — > u in C([0,T],V). Therefore 
u = v G L°°([Q,T],H m (T 3 ) n V"). From (2.10), it is easy to see that f t v G 
L°°([0, T],H m (T 3 ) n V), and thus by the Sobolev embedding theorem, we have u G 
C([0,T],F m (T 3 )nV). Again by (2.10), we see that u G C 1 ([0,T],i7 m (T 3 )ny). □ 

4. A Gevrey Class Regularity (Spatial Analyticity) 

Next, we show that (1.4) (or equivalently (2.10)) has additional regularity, un- 
der the assumption of the relevant regularity of u m . Namely we show that it has 
a specific type of Gevrey regularity which is analytic in space. In accordance with 
the discussion in Section 3, we do not expect solutions for initial data which is 
not analytic to become analytic instantaneously, unlike in the case of parabolic 
equations. (However, see [29], which proves the analytic regularity of the attractor 
of the Navicr- Stokes- Voigt equations, i.e., in the presence of viscosity and an ana- 
lytic forcing term. See also [49] for similar results concerning second-grade fluids.) 
The concept of Gevrey regularity was first used in the context of the Navicr- Stokes 
equations in [26] and was expanded to more general non-linear parabolic equations 
in [23] (see also [9,10]). For the three-dimensional Navier- Stokes- Voigt equations, 
Gevrey regularity for the attractor was shown in [29]. Gevrey regularity for the 
Euler equations in bounded domians has been studied in [4] and [34]. 

We define the Gevrey class of order s > for a given r > as the domain of a 
particular class of operators parameterized by r > 0: 

V{A r e TAlH23) ) := [u G H r (T 3 ) n V : Ue^^u^^) < oo} . 

For our purposes, we will work in the case s = 1, where T)(e rA 1 ) corresponds 
to the set of real analytic functions C"(T 3 ). For more on Gevrey classes, see 
[23,25,26,34,38,51]. We will need the following lemma, contained in [38] (see 
also [44]). 

Lemma 4.1. For ip,ip defined on T™ such that %jj G T>{A r / 2+l ^e rA 1 ) and tp G 
V{A r / 2+1 l 2 e TAl ' 2 ), where r>^andV-^ = 0,we have 




where c±, C2 > depend only on r. 

We are now ready to state and prove the main theorem of this section. The proof 
given in the present work is established rigorously using the Galerkin method. To 
establish our estimates, we follow some of the ideas of [38], which were applied 
formally to the 2D lake equations to establish formal a priori estimates. However, 
the results of [38] can be justified rigorously by using methods similar to those 
employed in the proof below. 
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Theorem 4.2. Let u in e V{A { - r+1 ^ 2 e' yAl/2 ) for some r > 3, r G N 7 and a > 0. 
Then there exists a unique solution u of (2.10) (with f = 0, for simplicity) satisfying 

(4.1) a^A^e^ Al/ \u(t)\\l 2(J3) + \\A*'*e^\{t)\\h m 

< a 2 \\A r ' 2 e° Al '\u m \\l 2 „ 3) + \\A r ^e- Al/2 u m \\ 2 L2{T3) +2 Cl f* k 3 (£)C 

Jo 

for all t G (—00,00), where r(t) := crexp^— C2 /j*' h(£) d£j , h > is defined by 
(4.4) below, c\,C2 > are given by Lemma 4-1 above, and «;(£) := C Q (1 + |£|) p ^, 
where p(r) = 5 (f) r — 1- 

The above theorem shows that for all t G (—00,00), u(-,t) is bounded in the 
space T>{A^ r+1 ^ 2 e T ^ A 1 ). Therefore, u is analytic with respect to the spatial 
variable. Furthermore, the radius of analyticity of the solution, which is bounded 
from below by r(i), may be shrinking as time increases, but it never collapses to 
zero in finite time, since r(£) > for all t. 

Proof. First, we observe that the case a = is already covered by Theorem 3.1; 
therefore, we assume that a > 0. Choose T > arbitrarily. Here again, we focus on 
the interval [0,T]. The proof for [-T, 0] is the same. Clearly V{A^ +1 ^ 2 e aAl ' 2 ) C 
£T r+1 (T 3 ), so the conclusions of Theorem 3.1 hold. It remains to prove (4.1), 
which we establish using the Galerkin approximation method. Consider the system 
(3.15). As discussed in Proof 2 of Theorem 3.1, this system has a solution un in 
C 1 ((— 00, 00), Hn) where Hn := span {fk}^ = i and (p^ are defined in Section 2. We 
denote ^ by f. Applying A T l 2 e TA 1 to both sides of (3.15a) and using the identity 

*l/2 d d tA 1/2 . 1 /2 tA 1/2 

e — UA r = — (e rA u N )-TA 1 ' 2 e TA u N , 

we have 



n 2 



+ ±(A r / 2 e TAl/2 u N )^fA^ +1 ^ 2 e TAl/2 u N = -A r ' 2 e TAl/2 P N B{u N ,u N ). 
dt 

Taking the L 2 inner product of this equation with A r / 2 e TAl/2 un and using the facts 
that A is a positive self-adjoint operator and V • un = gives 

~ (a 2 \\A^ Al/2 Vu N \\ 2 L2{T3) + \\A^ Al/2 u N \\ 2 L2{T3) 

(4-2) = a 2 T\\A r/2+1, ^ Al '\u N \\ 2 L2(J3) + T\\Ar/ 2 ^/^ Al ' 2 u N \\ 2 L2(J3) 

- (A r / 2 e TAll2 B(u N ,u N ),A r ' 2 e TAl ' 2 u N ) . 

We now estimate the last term in (4.2). We use techniques similar to [38,44]. Let 
us set tp = tp = un in Lemma 4.1 so that we have 



[A r ' 2 e TAU2 B{u N , u N ), A r l 2 e TAl/2 u N ) 



< Cl \\A r/2 UN\\l. +rc 2 \\A r / 2 e- Al/2 u N \\ L2 \\A r / 2+1 ^ Al/2 u N \\ 2 
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Using this in (4.2) gives 

— (a 2 ||^/ 2 e^ 1/2 V^||| 2(T3) + \\A^ Al/2 u N \\l 2{ri 
(4.3) < Cl ||^/ 2 UJV ||3 2(T3) +ra"\\A r / 2 +^ Al/2 Vu N \\l 2{T3) 

+ (i + rc 2 \\A r /\^ Al/2 u N \\ L . m ) \\A r / 2+1 ^ Al/2 UN \\l 2{T3) . 

Note that from the proof of Theorem 3.1, we have || J 4 r / 2 itjv(<:)||£2( T 3- ) < C a (l + 
\t\)P( r ) ■= K (t), where p is defined in the statement of the theorem. Recall that we 
defined r(t) := crexp ( -c 2 J * dA . Note that r(0) = a. The function h(t) > 



is defined by 

(h(t)) 2 :=l + a 2 ||^/V Al/2 V W m |l! 2( T3) 
(4.4) ,|t| 

+ \\A r ^e aAl/2 u m \\ 2 L2{T3) + 2 Cl / K 3 {8) M. 



This choice of r follows ideas from [38,44]. Now, since r, h > on [0, oo), we have 
r(t) = —C2h(t)T(t) < 0. Furthermore, we claim that 

(4.5) ||^/ 2 e T « Al/2 Ujv ||^(T3) < fc(t) 

for all time t > 0. Clearly (4.5) holds at i = 0, and therefore by continuity it holds 
for a short time. Let t* := sup{6> > such that (4.5) holds on [0,6)}. If t* = oo, 
we are done. Thus, suppose that t* < oo. Using the above bounds, integrating 
(4.3) on [0, t*], and recalling that r(0) = er, we have 

< a 2 \\A r / 2 e° Al/2 \7P N u m \\ 2 L2m + \\ A^e*^ P N u m || 2 2(T3) + 2 Cl f K 3 {9)d6 

Jo 

< h 2 (t*). 

Therefore the (strict) inequality (4.5) holds on [0,f*], and can thus be extended in 
time beyond t* , contradicting the definition of t* . Hence, the assumption that t* < 
oo must be false. Thus for all time, f + TC2\\A r l 2 e TAl ' 2 Mjv||z, 2 (T 3 ) < r + rc-xh = 0. 
Combining the above observations with (4.3) yields 

(4-6) — (a 2 ||^/ 2 e^ 1/2 V^||l 2(T 3 ) + H^V^ V|| 2 2(T3) ) < c lK 3 (t). 

Integrating (4.6) and using the fact that r(0) = a, we have 
(4.7) 

a 2 \\A r / 2 e^ Al/2 yu N (t)\\ 2 L2{T3) + \\A r / 2 e^ Al/2 u N (t)\\ 2 L2{T3) 

< a 2 \\A r / 2 e° Al/2 VP N u m \\ 2 L2{T3) + \\A r / 2 e^ 1/2 P N u m \\ 2 L2{T3) + 2c, f K 3 (£)^. 



o 



Set R := a 2 \\A r / 2 e° Al/2 Vu m \\ 2 L2{T3} + \\A r / 2 e° AU2 u m \\ L2{T3) and fix T e (0,T] 
arbitrarily. Define T m - m := min 0<t< j;r(t) = t(T) (since r is decreasing). Notice 
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that T m in depends on r, a, T, and R alone. Consider the Banach space B a ~ 
(which is in fact a Hilbert space) given by 

B« a ~ R :={^EH: a^A^e^^V^ 2 ^ + ||^/ V^%||£ 2(T3) < ^} , 
endowed with the indicated norm. Due to (4.7) and the fact that t(0) = a, we have 
a 2 \\A r / 2 e^ Al/ \u N (t)\\l 2{Ti} + \\A r / 2 e^ Al/2 u N (t)\\l 2m 

< sup (a^A^e^ Al/2 yu N (t)\\l 2{T3) + \\A^e^ AV \ N {t) ||| 2(T3) ) 

0<t<T 

< a 2 \\A r l 2 e° Al '\p N u m \\ 2 L2 „ 3) + \\A r / 2 e° Al/2 p N u m \\ 2 +2 Cl [ T 

Jo 

< a 2 \\A r / 2 e° Al/2 Vu in \\ 2 L2 „ 3) + \\A r / 2 e° Al/2 U in \\ 2 L2 „ 3) + 2c x C k s 

Jo 

Thus, ujv is uniformly bounded with respect to N in C([0,T],B Q We now 

find a uniform bound for the time derivative. To do this, we use the fact that 
||itij||_Ba _ <C , q ,||u||b° _ IMIb" - for s > 1/2, which was shown in Lemma 1 

s,tr,T,R s.cr,T,R. s.cr.T.R 

of [23]. Applying (I + a 2 A)' 1 to (3.15a), we calculate 
d 



dt UN 



= + a 2 A)- 1 P N B(u N ,u N )\\ Ba 

ga _ r+l,<r,T,R 



< C a \\B(u N ,u N )\\ B a _ <C a \\u N \\ Ba _ IIVujvIIj 



\B<* ^ ||".iV \\B" II v \\ B ° 

r-l,a,T,R r—1, a,T,R r 

|2 



< C a \\u N \\ Ba _ ||ujv|| B c _ < C a \\un\\b* 

r-l,(T,T,R r,v,T,R r,(T,T,R 

Thus, due to the uniform boundedness of the sequence ||itAr|| B a _ , we have that 

t .{j .T , R 

^ujv is uniformly bounded in C([0, T], a f R ) with respect to N, so that 
{"jv}j/=i i s equicontinuous with values in B a ~ . Let 0<f<l. Thanks to the 
compact embedding B a ~ <-} B a ~ and the Arzcla-Ascoli Theorem, we see 

r,cr,T,i? r — e,cr,T,R 

by the above bounds that there exists a subsequence of {wjv}^ =1 , which converges 
in C(\0, T] , B a ~ ) to an element w £ C(\0, f],B a ~ ). As shown in Proof 

2 of Theorem 3.1, this sequence also converges in C([0, T], V) to the unique solution 
m given in Theorem 3.1. By the uniqueness of limits, u = w £ C([0, T], -B" -^J, 
and furthermore, we have the bound 

sup \\u{t)\\ 2 Ba _ 
te[o,T] r " e ' '■ T - R 

< a 2 \\A r / 2 e° Al/2 Vu m \\ 2 L2{T3) + \\A r ' 2 e° Al/ \*»\\l 2(T 3) + 2 Cl [ T 

Jq 

Since this bound holds uniformly for all < e <C 1, it also holds for e = 0, which 
can be seen by using the Fourier series representation 

!K*)lll« _ = £(i + a 2 |j| 2 )|j| 2 ^>e 2 -- ^uM 2 , 

r-e,a,T,R '—^ 
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taking the limit as e — > and passing the limit inside the infinite sum using, for 
example, the Lebesgue Monotone Convergence Theorem for the counting measure. 
Here, we have used the notation Uj(t) for the Fourier coefficients of u{t). Thus, 

a 2||^/2 e WV 2vu(t)|| 2 2(T3) + || ^r/2 e W^\ W ||2 2(T3) 

< a 2 ||^/ 2 e ^ 1/2 V U »||| 2(T3) + P'VV^V-Hi +2ci f T ^(Odt 

Jo 

for all t £ [0, T]. Finally, by choosing t = T in the above inequality, and by recalling 
that r min = t(T) and that T was chosen arbitrarily in (0, T], we have established 
(4.1) for u, with T playing the role of t. □ 

5. Convergence to the 3D Euler Equations and a Blow-up Criterion 

We now consider the three-dimensional Euler equations, namely (1.1) with v = 
and periodic boundary conditions, or equivalently (2.9). It is known that if u m £ 
H S (T 3 ) nV for s > 5/2, then there exists a unique solution to these equations 
in C([0,T*],iF(T 3 )nF)) nC 1 ([0,T,],P- 1 (T 3 ) n V)) for some T» > (see, e.g., 
[31,32,39,42,43,53]). Observe, however, that the question of existence of weak 
solutions for the Euler equations, i.e., for u m £ H s , s < 5/2, or u m £ C s , S £ [0, 1] is 
still open. In fact, the only results for short-time existence are those for u m £ C 1 ' 6 , 
proven in [39], or u m £ H s , s > 5/2, (see, e.g., [31,42,43]), where the solutions 
are well-posed in the sense of Hadamard. Recently, it was shown in [5] that the 
three-dimensional Euler equations are ill-posed in the space C s . 

For arbitrary u m £ H s (T n ) n V, with s > 3, we consider the maximal time 
interval [0, T max ), T max > T*, for which a solution to the Euler equations exists and 
is unique. It is a major open problem in mathematics to determine whether T max 
is finite; that is, whether solutions exist globally in time or experience blow-up in 
finite time (see the recent surveys [6] and [16]). In the present section, we present a 
criterion for blow-up of the Euler equations using methods similar to those in [33] . 
There are few other known criteria for blow-up of the Euler equations which are 
directly checkable for the Euler equations themselves (see, e.g., [7,17,19,22]). 

Theorem 5.1 (Convergence). Given initial data u m , u™ £ H S (T 3 ) for some s > 3, 
let u,u a £ C7([0,T],iJ s (T 3 )) n C 1 ([0, T], i/ s_1 (T 3 )) be the corresponding solutions 
to (1.1) and (1-4), respectively, where < T < T max , andT max is the maximal time 
for which a solution to the Euler equations exists and is unique. Then: 

(i) For all t £ [0,T], 

- u a (t)\\ 2 L2{T3} + a 2 \\V(u(t) - u a (i))||£2(T3) 

(5.1) < (IK 1 - <||! 2(T3) + a 2 ||VK" - <)||| 2(T3) )e ct + Ca 2 (e ct 1). 
(ii) Consequently, if 

(5.2) |K" - <||! 2(T 3) + a 2 ||V( U ™ - 0||| 2(T3) as a -+ 

(in particular, when it™ = u m for all a > 0), then u a — >• u in C([0,T], H) , 
as a —> 0. 

Remark 1. Regarding part (ii) of Theorem 5.1, one may find many sufficient 
conditions for (5.2) to hold. For example, if u™ — > u m in L 2 (T 3 ), as a — > 0, and 
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either a||w™||#-i( T 3) — ► 0, as a — ► 0, or 1 1 it™ 11^(13) < C/a 17 for some a £ (l,s], 
then (5.2) holds, due to the interpolation inequality (see, e.g., [1]) 

ii in in\\ ^ s~t\\ in mill — 1/cr m in in\\l/& 

IK ~ u Whht 3 ) < C\\u a -u \\ L *fa\\u a - u ||^ (T 3)- 
Therefore u Q — >• u in C([0,T],_ff), as a -> 0. 

Remark 2. With slightly more effort, one can prove similar convergence results 
to Theorem 5.1 for it m ,u™ £ H S (T 3 ) with s > 5/2. However, for simplicity of 
presentation, we take s > 3. In fact, in lieu of the result in [48], which shows the 
local well-posedness of the Euler equations in the Besov space i(T"), n > 2, one 
can prove a similar result to Theorem 5.1 assuming only that u m , u™ £ B^ i(T"). 

Proof. Subtracting (1.4) from the Euler equations gives 

d t {u - u a ) + a 2 d t Au a + V(p - p a ) 

= — (u ■ V)lt + (u a ■ V)u a 

= — ((it — U a ) ■ V)ll + ((u — U a ) ■ V)(lt — U a ) — (u • V)(li — U a ). 

We note that since u,u a £ C' 1 ([0, T], H 11 " 1 ) and s > 3, we have itt, <9tit a , Adtu, 
AdtUa £ C([0,T], L 2 (T n )), so that taking inner products with these terms is jus- 
tified. This is because i? s_1 (T 3 ) is an algebra for s > 3 (see, e.g., [1]), so we 
also have u ■ Vu, (u a ■ V)u a £ C([0, T],H s ~ 1 (T n )). Using (2.6) and the fact that 
V • u = V • u a = 0, we integrate against u — u a to find 

n~r\\ u ^ u a\\ 2 L2(T3\+ a 2 [ AdfU a • (u — u a ) dx = — [ ((u — u a ) -V)u ■ (u~ u a ) dx. 
2 at J T 3 J T 3 

The exchange of the order of the integral and the time derivative is justified here 
and below since it, u a are C 1 in time with values in H s . Adding and subtracting 
the term a 2 | T3 (A« ( ) • (it — u a ) dx, we have 

ollll" ~ u a\\h(T3) - ® 2 / (Au t - Ad t u a ) ■ (u - u a ) dx 
i at j T 3 

,2 



+ a / (Aui) ■ (u — u a ) dx = — / ((it — u a ) ■ V)it • (it — u a ) dx. 
Jt 3 Jt 3 

Integrating by parts and rearranging gives 

(5.3) = — a 2 / (Au t )(u — u a ) dx — / ((u — u a ) - V)u • (u — u a ) dx. 



it 3 Jt 3 
Using the fact that it satisfies the Euler equations, we have 

—a 2 / (Ait t ) ■ (it — u a ) dx — a 2 I (—ut)-A(u — u a )dx 
Jt 3 Jt 3 

— a 2 / ((it • V)it) • A(it — u a ) dx 
Jt 3 



a 2 



-a 2 



/ Vp-A(u — u a )dx 
Jt 3 

V((it • V)it) • V(it - it Q )dx. 
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The last equality follows by using integration by parts on the first integral and 
noting that the second integral term is zero because of (2.8) and the fact that 
div(A(u — u a )) = 0. Thus the right-hand side of (5.3) becomes 

—a 2 / (V((tt ■ V)ii))V(u — u a ) dx + / ((it — u a ) ■ V)it • (it — u a ) dx. 

Since we are assuming u is a regular solution, i.e., u G H s for some s > 3, the H 3 
norm is finite, so the second term is bounded above by 

||VM|| L =C(T3)||U-U Q ||^2 (T 3) < C||lt|| ff 3( T 3)||u - U Q ||i2( T 3) < C\\ U - U a || i 3 (T3) ■ 

Using the Cauchy-Schwarz inequality on the first term gives 
-a 2 / (V((u ■ V)u)) • V(« - u ) *e 

JT 3 

2 



< q ||V((u • V)u)|| L 2 (T 3 ) ||V(m - it Q )|| L 2(T 3 ) 

< a 2 (|||Vu| 2 |U2 (T3) + p 2 U || L2(T 3)) ||V(« - U a )\\ L 2 

.2 



C2U 2 . 



< a 2 (||VU||L4( T 3) + ||u||iJ2( T 3)) ||V(li - U a )\\ L 2 

< a 2 (C||Vu||#2( T 3) + ||«||h2(t3)) ||V(u - U a )\\ L 2 
<2Ca 2 \\V(u-u a )\\ L2 
<C 2 a 2 +a 2 \\V(u-u a )\\ 2 L 2. 

Collecting the above estimates we have from (5.3) 

4 (\\ u - u a\\h(Ti) +" 2 ||V(u-u. Q )|| 2 j2(T 3 ) ) 

(5.4) <d (||«-U a ||i aCl » ) +Q! a ||V(«-« a )||i 3fia) ) 

Gronwall's inequality then gives 

\\u(t) - u a (t)f L2(J3) + a 2 \\V(u(t) - u a (t))||| 3(x3) 

(5.5) < (\\u™ <|j 2 2(T 3) + a 2 \\V(u m <)|| 2 2(T 3 ) )e Clt + ^{e c ^ - 1). 

This proves (5.1). Next, if <™ w m in L 2 (T 3 ) and a 2 || V<"|j 2 2(T3) < M < 00 
is bounded as a — > 0, noting that the left-hand side of (5.1) is bounded below by 
\\u — u a ||f,2(-T3)! we take the limsup a _,. to find that u a — >• it in C([0, T), iT). □ 

Due to Theorem 5.1, we have the following result. 

Theorem 5.2 (Blow-up Criterion). Assume u m £ H s , for some s > 3. Suppose 
there exists a finite time T** > such that the solutions u a of (2.10) with it™ = u ln 
for each a > satisfy 

sup limsupa 2 ||Vit a (t)j| 2 2( - T 3') > 0. 

te[0.T") a->0 

TTien the Euler equations with initial data u m develop a singularity in the interval 
[0,T**\. 

Proof. Suppose to the contrary that the solution it to the Euler equations is regular 
on [0,T**]. By (3.3) we have the "modified energy" equality 

(5.6) \\u a (t)\\ 2 L2(j3) + a 2 \\Vu a (t)\\h { ^ } = IklW 3 ) + a 2 \\Vu in \\l 2{r3) , 
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for all t > 0. Taking the limsup of (5.6) gives 

(5.7) \Ht)\\ 2 L2{T3) +limsupa 2 ||V UQ (t)||| 2(T3) = K n ||£ 2(T3) , 

a—>0 

for all t e [0,T**], since \\u a (t)\\ L 2 -> \\u(t)\\ L 2 for t G [0,T**] by (5.5). Since 
regular solutions of the Euler equations conserve energy, we have the energy equality 
||u(£)|| L 2( T 3) = ||u m || L 2 (T 3) for t £ [0,T**]. Thus (5.7) contradicts the hypothesis 
and u must blow-up in [0,T**]. □ 



6. The MHD-Voigt Case 

In this section, we consider a Voigt-type regularization of the inviscid, irressitive 
magncto-hydrodynamic (MHD) equations, given by 

(6.1a) ~a 2 d t Au + d t u + (u ■ V)u + Vp + ^V|B| 2 = (B ■ V)£, 

(6.1b) -a 2 M d t AB + d t B + (u ■ V)fl - (B ■ V)u + Vq = 0, 

(6.1c) V-i3 = V-u = 0, 

(6.1d) B(Q) = B m , u(0) = u"\ 

where a, «m > 0, and the boundary conditions are taken to be periodic, and we 
also assume as before that 

/ udx = I B dx = 0. 

Jt 3 Jt 3 

Here, the unknowns are the fluid velocity field u(x,t), the fluid pressure p{x,t), 
the magnetic pressure q{x 1 t), and the magnetic field B(x,t). (By formally taking 
the divergence of (6.1b), we find that V(? = 0, but this is not assumed a priori.) 
Note that when a = qlm = 0, we formally retrieve the inviscid, irressitive MHD 
equations, i.e., the case where the fluid viscosity and the magnetic diffusivity are 
equal to zero. The viscous Bardina model for the MHD case has been studied 
in [35], and another viscous MHD a- model has been studied in [40]. Here, we study 
the inviscid case, and hence prove stronger results than those reported in [35]. 

We first prove that this system has a unique local (in time) solution. Let us 
show that (6.1) has a unique short-time solution. To do this, we follow [11] and use 
the Picard-Lindelof Theorem. Applying P a (see Section 2) to (6.1), we obtain the 
equivalent system 

Q 

(6.2a) — (a 2 Au + u) = B(B,B) - B(u,u), 

at 

(6.2b) ^(a 2 M AB + B) = B(B,u) - B(u,B), 

(6.2c) B(0) =B m , u(0) =u m . 

Note that it is possible to recover p and q by using (2.8) (see, e.g., [41]). Denote 
v = (a 2 A + I)u, Z = (a 2 j A + I)B, N^v, Z) = B(B, B) - B(u, u), and N 2 (v, Z) = 
B(B,u) — B(u,B). Then 6.2 is equivalent to the system 

d (v\ {Ni(v,Z)\ fv(0)\ fv in \ f-a 2 Au in + u v 



^ 6 ' 3 ^ dt\z) \N 2 (v,Z)J' \Z{0)J \Z m )' \-a 2 M AB in + B i 
Using this form of (6.1), we will now prove the following theorem. 
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Theorem 6.1 (Short-time existence and uniqueness). Let v m ,Z m e V . 

Then there exists a time T = T(\\v m \\y , ||Z m ||y) > such that (6.3) has a unique 

solution (v,Z) G C l ([-T,T\,V'), or equivalently, (u,B) € C 1 Q-T,T],V). 

Proof. Wc follow almost word by word the last section of [11]. To show local 
existence, it is enough to show that N± and N 2 are locally Lipschitz in the space 
V'. By the Poincare inequality (2.2) and (6.2a) we have 

||i\Ti(i;i,Zi)-.Ni(«a J Z2)||v' 
= \\B(Bi,Bi) - - B(B 2 , ft,) + B(u 2 ,u 2 )\\ v , 

= \\B(ui -u 2 ,u 2 )-B(u 1 ,u 1 -u 2 ) + B{Bi -B 2 ,B 2 ) - B(B U Bi - B 2 )\\v 

< ||B(tii-« 2 ,ti2)-B(ui,tii-« 2 )||v' + ||-B(Bi-B2,B2)-B(Bi,iBi-ft)||v' 
= sup I (B(ui — u 2 ,u 2 ) — B(ui, u\ — u 2 ), w) I 

{weV:\\ W \\ = l} 

+ sup \(B{B 1 -B 2 ,B 2 )-B{B 1 ,B 1 -B 2 ),w)\ 

{w£V:\\w\\ = l} 

< c\ ui - U2 i i/2 k - " 2 ii i/2 imi + cki i/2 im 1/2 ik - U2 || 

+ C\B X - B 2 | 1 /2j| 6l _ B 2 fl 2 \\B 2 \\ + CiBil^HBill^HB! - B 2 || 

< CA- 1/4 || Ul - ua|| (|| Ul || + IKH) + CX- 1 '^ B 2 || (||Bi|| + ||B 2 ||) 
<2Ci?A- 1/4 (|| Ul - U2 || + ||B 1 -B 2 ||) 

< Ci?Ar 1/4 (||wi - «2||v + ||Zi - Z 2 \\v>) , 

where R is chosen so that ||iti||, ||u 2 ||, ||Bi||, ||B 2 || < R. Next, for N 2 , we have by 
the Poincare inequality (2.2) and (6.2a), 

WN^v^Zx) - N 2 {v 2 ,Z 2 )\\ v , 
= ||B(Bi,ui) - +B{u 2 ,B 2 ) + B{B 2 ,u 2 )\\ v , 

< ||fl(«i,Bi - B 2 ) - S( Ul - U2,B2)||v + ||B(Bi.ui - ua) - B{B\ - B 2 ,u 2 )\\ v , 
= sup I (B(ui,Bi -Ba) — B(ui —u 2 ,B 2 ),w) \ 

{w€V:\\w\\=l} 

+ sup \{B{B 1 ,u 1 -u 2 )-B{B 1 -B 2 ,u 2 ),w)\ 

{weV:\\w\\ = l} 

< cm^ihi^HBi - B 2 ]i + c|«i - u 2 | 1/2 K - « 2 || 1/2 ||b 2 || 

+ CIBxI^UBxH 1 / 2 ^! - Ua || + C\Bi - B 2 | 1/2 ||Bi - B 2 || 1,/2 ||u 2 | 

< CA- 1/4 (|| Ul || + HualDIIBx - B 2 || + CA^d^H + ||B 2 ||)H - u 2 \\ 
<2Ci?A- 1/4 (|| Ul - U2 || + ||B 1 -B 2 ||) 

< CR\; lfi (||«! - v 2 \\ v , + \\Z X - Z 2 \\ v ,) , 

where again R is chosen so that ||mi||, ||w 2 ||, ||Bi||, j|B 2 || < R. Thus, N\ and N 2 
are locally Lipschitz in V', so the right-hand side of (6.3) is as well. Therefore by 
the Picard-Lindelof Theorem, there exists a unique solution (v, Z) to (6.3) such 
that v, Z € C l {[— T, T],V) for some T > which may depend upon the initial 
conditions. □ 

We next show that in fact, we have global existence. 
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Theorem 6.2 (Global existence and uniqueness). Let v m ,Z m e V (equivalently 
u m ,B m £V). Then (6.3) has a unique solution (v,Z) £ C ((— oo, oo),V') (equiv- 
alently (u, B) G C 1 ((-oo, oo), V) ). 

Proof. It is sulBcient to show that on the maximal time interval of existence 
and \\Z(t) \\v are finite. Let [0, T*) be the maximal interval of existence. If T» = oo, 
we are done, so suppose T* < oo. Then we must have limsup t _;. Tt IkWIIv = 00 
or limsup t _;. Ti ||Z(i)||y = oo, otherwise we could use Theorem 6.1 to extend the 
solution further in time, contradicting the definition of T* . Therefore, we must have 

(6.4) limsup ||u(t)|| = oo or limsup = oo. 

t-¥T~ i->T~ 

Taking the inner product of (6.2a) with u and (6.2b) with B and integrating by 
parts, we have 

(6.5a) ^(a 2 \\u\\ 2 + \u\ 2 )=(B(B,B),u), 

(6.5b) ^-(a 2 M \\B\\ 2 + \B\ 2 ) = (B(B,u),B) = -(B(B,B),u), 

where in the last equation, we used (2.5). Adding (6.5a) and (6.5b) and integrating 
in time, we obtain on [0,T*) 

(6.6) a 2 \\u\\ 2 + a 2 M \\B\\ 2 + M 2 + \B\ 2 = a 2 \\u in \\ 2 + a 2 M \\B in \\ 2 + \u m \ 2 + \B m \ 2 . 

Since the right-hand side is finite for all time, this contradicts (6.4). The proof is 
nearly identical for (— T*,0]. Thus, T* = oo. □ 

The higher-order regularity for system (6.1) holds as well. 

Theorem 6.3. Let u"\B ln £ H m (T 3 ) n V, for m > 1. Then the unique solution 
(u, B) of (6.1) given by Theorem 6.2 lies in C((— oo, oo), H m PI V). 

Proof. The proof follows nearly the same steps as in Section 3, mutatis mutandis, 
so we only sketch the main ideas. The primary difference lies in handling the 
additional variables given by the magnetic terms. We mention that the work here 
is done formally, but can be made rigorous by following similar ideas to those in 
either of the two proofs of Theorem 3.1 given above. One must work inductively, 
first estimating the H m norms for m = 0,1,2,3 as in steps 1-3 of Proof 1 of 
Theorem 3.1, and then obtain the estimates for general to > 4, as in step 4. The 
cases to = 0, 1 are already given by (6.6). Here, we only show the to = 2 case, 
where the central issue is more transparent. The other cases are more complicated 
notationally, but not conceptually. We first apply to (6.1a) for an arbitrary j3 
with \P\ = 1, and integrate the result against d^u to obtain 

i£ (a 2 ||d^|j 2 + \d^u{ 2 ) = (B(d^B,B),d u) + (B(B,d^B),d^u) 

- (B(d p u,u),d^u) - {B(u,d e u),d^u) 

(6.7) = (B(d B, B), d u) + (B(B, 8^B),d^u) 

- (B(d u,u),d^u) 
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Where we have used (2.6). Next, we apply d 13 to (6.1b) and integrate the result 
against d^B 

~ (a 2 M \\d^B\\ 2 + \d?B\ 2 ) = (B(d^B,u),d^B) + (B(B,d^u),d^B) 

- (B{d^u,B),d^B) - (B(u,d B),d^B) 

(6.8) = (B(d B, u),d p B) + (B(B, d u), d B) 

- {B(d^u,B),d^B) 

Since (B(B,d^u),d^B) = -~(B(B,d^B),d^u) by (2.5), we may add (6.7) and (6.8), 
to obtain a very important cancellation of the terms that involve the highest order 
derivatives. This gives us 

\± (a 2 M \\d^B\\ 2 + \d?B\ 2 + a 2 \\d?u\\ 2 + \d^u\ 2 ) 

(6.9) = (B(d B, u), d^B) - (B(d u, B), d^B) 

+ (B(d B, B), d u) - (B(d u, u), d^u) 
^SUVBII^^IIVSIUa^llVulUe^) 

+ l|V-u|| L 2 (T 3 ) ||Vu|| L 3 (T 3 ) ||Vu|j i 6 (T 3 ) 

< C\\B\\ ||B|| ff 3/2 (T3) ||u||h2(t3) + C\\u\\ \\u\\ H3/ 2 {T3) \\u\\ h »(t*) 

< C a ,a M (||S||#3/2( T 3) + II M II^3/2(t3) + IMIff2(T 3 )J 

since ||£>||, ||u|| are uniformly bounded from (6.6). Summing over all /3, with \[3\ = 1, 
and using the norm equivalence (2.3) yields 

\j t {o? u \\Bf HHT3) + \\B\\ 2 + a 2 ||u|| 2 H2(T3) + \\u\\ 2 ) 

(6.10) < C a>aM (\\B\\ 2 H . m + \\u\\ 2 H2{ri) ) . 

Thus, the H 2 (T 3 ) norms have at most algebraic growth by Gronwall's inequality. 
To complete the proof, the estimates for the higher derivatives are carried out in 
a similar way, using the methods of steps 3 and 4 of Theorem 3.1, and making 
use of the analogue of the important cancellation that was used to obtain (6.9). 
We emphasize again that the estimates here are derived formally, but can be made 
rigorous, as described above. □ 

Remark 3. One can also obtain results for equations (6.1) concerning the Gevrey 
regularity, similar to those obtained in Section 4 for equations (2.10). However, 
these will be omitted here because the concepts are exactly the same, and one only 
needs to combine the ideas of the proofs of Theorems 4.2 and 6.3 to recover them. 
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